Quantum quench dynamics of the sine-Gordon model in some solvable limits 
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In connection with the the thermalization problem in isolated quantum systems, we investigate 
the dynamics following a quantum quench of the sine-Gordon model in the Luther-Emery and 
the semiclassical limits. We consider the quench from the gapped to the gapless phase as well as 
reversed one. By obtaining analytic expressions for the one and two-point correlation functions of 
the order parameter operator at zero-temperature, the manifestations of integrability in the absence 
of thermalization in the sine-Gordon model are studied. It is thus shown that correlations in the 
long time regime after the quench are well described by a generalized Gibbs ensemble. We also 
consider the case where the system is initially in contact with a reservoir at finite temperature. The 
possible relevance of our results to current and future experiments with ultracold atomic systems is 
also critically considered. 

I. INTRODUCTION 

Recent experiments in the field of ultracold atomic gases have spurred much interest in understanding the ther- 
malization dynamics of isolated quantum systems. [1-28] So far, the latter were considered mere idealizations of real 
systems, as most of the many particle systems of interest to quantum statistical mechanics, such as solids and quan- 
tum fluids, are strongly coupled to their environments. However, the creation of large ensembles of ultracold atoms 
with highly controllable properties, which remain fully quantum coherent for relatively long times (compared to the 
typical duration of an experiment), has completely changed this perception. This has also raised concerns about the 
mechanisms of thermalization in these systems, and even in some recent experiment [6] lack of thermalization has 
been observed. 

The problem of thermalization in isolated quantum systems can be also posed as the study of the dynamics of the 
system following a quantum quench. That is to say, the study of the response of the system to a change of a control 
parameter of the Hamiltonian over a time-scale much shorter than any other relevant time scale of the system, so 
that the sudden approximation can be applied. Therefore, it is assumed that for t < the Hamiltonian is Hi and the 
system is in a given eigenstate of it, |<f>;). At t — the Hamiltonian is changed to Hf and thus for t > the system 
evolves unitarily in isolation according to the dynamics dictated by Hf. The quantum quench can be also used to 
describe the evolution of a system that has been prepared in a given state that is not an eigenstate of the Hamiltonian. 
Thus, the question naturally arises is whether a set of sufficiently interesting obscrvablcs of the system reach some 
form of stationary state that can be described by a standard Gibbs ensemble. In such a case, we would speak of 
thermalization to a standard statistical ensemble (microcanonical, canonical, or grand-canonical). However, this may 
not be the case as it turns out that for several integrablc models [11, 14, 15, 19, 26] the long-time behavior of certain 
observables can instead obtained from a generalized Gibbs ensemble characterized by a different temperature for 
each eigenmode of the Hamiltonian Hi [11, 14, 26-28]. The non-standard or generalized Gibbs ensembles follow from 
maximizing the von Neumann entropy with the constraints imposed by the set integrals of motion larger than the total 
energy and particle number. As to the situation concerning non-integrable systems, the issue of the thermalization 
dynamics is still not fully understood although some recent results [20, 21, 24] indicate that thermalization to a 
standard ensemble should eventually occur at sufficiently long times. However, the actual thermalization dynamics 
and how it depends to how close the system is to integrability still remain very poorly understood. Indeed, as far 
as one dimensional systems are concerned (for which the strong kinematic constrains usually lead to integrability 
being more ubiquitous than in higher dimensions), numerical simulations have found conflicting results concerning 
the existence of thermalization [12, 13, 23]. Indeed, in a recent work, Rigol and coworkers [29, 30] has also pointed 
out that the statistics of the constituent particles may also play an important role in determining the thermalization 
dynamics. 

Furthermore, in addition to studying the quench dynamics starting from a pure state, it can be also interesting to 
consider quantum quenches where the system is initially prepared in a thermal state by being at t < in contact 
with an energy reservoir at temperature T. After the quench, the system is isolated from the reservoir and evolves 
unitarily according to the dynamics dictated by different Hamiltonian Hf. This is the simplest kind of mixed initial 
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state one can consider and it allows to analyze the effect on the quench dynamics of mixing (by means of the initial 
temperature) a fraction of the excited states with the ground state of . 

In this work, we shall analyze the quench dynamics of the sine-Gordon model (sGM). This is an intcgrable field 
theory, but its exact solution is indeed quite difficult to deal with as the elementary excitations satisfy a rather 
non-standard algebra. Instead, we confine ourselves to two limits in which the model can written as a quadratic 
Hamiltonian. In one of these two limits, the so called Luthcr-Emcry (LE) limit, using a trick called refermioniza- 
tion, [31, 32] the model describes a system of one-dimensional massive (i.e. gapped) Dirac fcrmions. The fermions and 
their anti-particles (or holes, to use a solid-state Physics language) describe the solitonic and anti-solitonic excitations 
of the sGM, which, in the LE limit, happen to be non-interacting. In another limit, the model can be rather well 
approximated by a quadratic model of massive bosons. The latter describe a series of bound states of solitons and 
anti-solitons in the limit where interaction between them is strongly attractive. As explained elsewhere, [11, 26] the 
quench dynamics of a quadratic Hamiltonian can be solved exactly via a time-dependent Bogoliubov transformation. 
We shall consider here two kinds of quenches, which correspond to the appearance and disappearance of the mass 
term (i.e. the gap) in the sGM. At zero temperature, we found correlations when the system is quenched from the 
gapped to the gapless phase at zero or low temperature exhibit an exponential decay with a correlation length/time 
fixed by the gap. This result is in agreement with the results of Calabrese and Cardy, [9, 15] based on a mapping to a 
boundary conformal field theory At high temperatures, however, the correlation length (time) is determined by the 
temperature. At intermediate temperatures, the system will exhibit a crossover between the zero-temperature (gap 
dominated) and the high-temperature (temperature dominated) regimes. On the other hand, correlations following 
a quench from the gapped to the gapless phase at the Luther-Emery and the semi-classical limit exhibit a somewhat 
different behavior, which may indicate a break-down of the semiclassical approximation or a qualitative change in the 
behavior of correlations as one moves away from the Luther-Emery limit. 

This paper is organized as follows. In section II, we introduce the (quantum) sine-Gordon mode, its phases, and 
describe the problem of quantum quenches in this model. In section III, we consider quenches in the Luther-Emery 
limit at zero temperature, whereas in section IV we take up the quenches in semiclassical limit. In section V, we 
consider the effect of thermal fluctuations in the initial state assuming that contact with the energy reservoir is 
removed at the time when the system is quenched and therefore, it subsequently evolves in isolation. In section VI, 
we take up the issue of the long-time asymptotic behavior of correlations and expectation values, and show that it 
is given by a generalization of the Gibbs ensemble, as recently pointed out by Rigol and coworkers. [14] Finally in 
section VII we discuss some possible experimental consequences of this work. We provide a summary of the conclusions 
in section VIII. The results for the correlations at zero temperature that are reported in this manuscript have been 
available previously as a part of an unpublished preprint. [28] Since our preprint appeared other authors have also 
considered the thermalization dynamics in the sine-Gordon model. In particular (see this special issue), Sabio and 
Kchrcin [33] used a flow equation method. Fiorctto and Mussardo used [34] form factor methods to tackle quantum 
quenches in the sine-Gordon model, finding strong evidence that the long time behavior of a local operator is described 
by a generalized Gibbs ensemble. Furthermore, using the time-dependent renormalization group [35, 36], Barmettler 
and coworkers [37] have investigated the quench dynamics of the XXZ model, which in the continuum limit reduces 
to the sine-Gordon model 



where : . . . : stands for normal order of the operators, [31, 32, 38] ao is a short-distance cut-off, and the phase and 
density fields, 9(x) and 4>(x), are canonically conjugated in the sense that the obey [4>{x), d x '9(x')] = iirS(x — x'). 
This model can be regarded as a perturbation the Luttinger model (see e.g. Refs. 31, 32, 38 for a review), which 
still yields an integrable model. In equilibrium the model is known to exhibit two phases, which, according to the 
renormalization group analysis, [31, 32] and for infinitesimal and positive values of the coupling in front of the cosine 
term, roughly correspond to K < 2 (gapped phase) and K > 2 (gapless phase). 

In order to study the non-equilibrium (quench) dynamics, we will consider two different types of quenches: the 
quench from the gapless to the gapped phase and the reversed process, from the gapless to the gapped. In the first 
case, we assume that the dimensionless coupling g(t) is suddenly turned on, i.e. git) — g6(—t). With this choice, 



II. THE SINE-GORDON MODEL 



The sine-Gordon model is described by the following Hamiltonian: 




(1) 



(2) 
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Hi = H sG (t < 0) is a Hamiltonian whose ground state exhibits a frequency gap, m, to all excitations, whereas 
Hi = H sG (t > 0) has gapless excitations. Conversely, in the second case, we consider that g(t) is suddenly turned 
off, i.e. g(t) = g0(t). In this case, the ground state of Hi is gapless whereas the Hamiltonian performing the time 
evolution, Hf , has gapped excitations. However, although both Hi and Hf define integrable field theories, for a general 
choice of the parameters K and g, the quench dynamics cannot be analyzed, in general, by the elementary methods 
of Refs. [11, 26]. Nevertheless, in two limits, the Luther-Emery (which corresponds to K = 1, see Sect. Ill) and in 
the semiclassical (that is, for K <C 1, Sect. IV) limits, it is possible to study the quench dynamics by methods similar 
to those of Ref. 26. However, as explained above, the statistics of the elementary excitations happens to be different 
in these two limits. 



III. THE LUTHER-EMERY LIMIT 



A. Introductory remarks 



Let us start by considering the sine-Gordon model, Eq. (3), for K = 1, which is the so-called Luther-Emery limit. 
It is convenient to introduce rescaled density and phase fields, which will be denoted as (fi(x) = K~ l / 2 <j)(x) and 
(p(x) = K 1 ^ 2 9(x). Thus, the Hamiltonian of Eq. (1) becomes: 



H sG (t) = ^ I dx: (d x tp) 2 



m ~n2 hvg(t) 
(d x (p) : j- 



dx COS Kip, 



(3) 



where k = 2\JK. At the Luther-Emery limit k, = 2 (i.e. K = 1) and the model can be rewritten as a one- 
dimensional model of massive Dirac fermions with mass by using the following bosonization formula for the Fermi 
field operators, [31, 32] 



ip a (x) 



g a e 



is a TT /4 



«4>«(x) 



'2na 



(4) 



where s r = — si = +1 and the chiral fields <j> r (x) = <p(x) + <p(x), and <fii(x) = <fi(x) — <p(x). For computational 
convenience, we choose the Majorana fermions in Eq. (4) to be r\ r = a x and r)i = ia y , where a x and a y are the familiar 
Pauli matrices. In addition, we note that the gradient terms in Eq. (3) can be written as the kinetic energy of free 
massless Dirac fermions in one dimension: [31, 32, 38] 



ihv \ dx : ipl(x)d x i/j r (x) - ipj (x)d x ipi(x) 



H () 



As far as the cosine operator is concerned, the bosonization formula, Eq. (4), implies that: 

ipl(x)ipi(x) + ijjj (x)i/j r (x) = cos2ip(x) 

7ra 



r 

7T 



cos2ip(x) : 



(5) 

(6) 
(7) 



where T = ia x a y . This is almost the cosine term of the sGM in the LE limit (cf. Eq. 3), except for the presence of 
the operator T. However, we note that r 2 = 1 and that this operator also commutes with Hq and with operator in 
the left hand-side of Eq. (6). The first property implies that the eigenvalues of T are ±1 whereas the second property 

implies that H^it) = Hq + hvg(t) J dx tpl(x)tpi(x) + (x)ip r (x) and T can be diagonalized simultaneously. Upon 

choosing the cigenspace where T = — 1, we obtain that 

#le(*) = -itw / dx : ipl(x)d x ip r (x) - ip\ (x)d x ipi(x) : 



+ hvg(t) j dx [^(aMO) +^(20^0)] , (8) 

is equivalent to Eq. (3) when k = 2. 

To gain some insight into the phases described by the sGM, let us first consider the Luther-Emery Hamiltonian, 
i?LE m two (time-independent) situations: i) g(t) = (the gapless free fermion phase, which coincides with the 
Luttingcr model for K = 1 [31, 32, 38]), and ii) g(t) = g > 0, i.e. a time-independent constant (which corresponds to 
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the gapped phase). In order to diagonalize the Hamiltonian, it is convenient to work in Fourier space and write the 
fermion field operator as: 



(9) 



where a = r, I. The limit where the cut-off do — > + should be formally taken at the end of the calculations, but in 
some cases we shall not do it in order to regularize certain short-distance divergences of the quantum sGM. It will also 
useful to introduce a spinor whose components are the right and left moving fields, and which will make the notation 
more compact: 



*(p) 



Iprip) 
Iplip) 



, H(p) 



fkv (p) 
-hwo(p) 



Thus the Hamiltonian for the gapless phase, Ho, reads: 



(10) 



(11) 



where ujq(p) = vp is the fermion dispersion. However, the Hamiltonian of the gapped phase, corresponding to 
g(t) = g > 0, Hue,, is not diagonal in terms of the right and left moving Fermi fields. In the compact spinor notation 
it reads: 



H LE = j2-^(p)-np)-^(p)-, 



(12) 



where 



H{p) 



hujQ (p) hm 
hm -~tvjj(){p) 



being m = vg. Nevertheless, i?LE can be rendered diagonal by means of the following unitary transformation: 

Hp) 



i>c(p) 




i>v{p)_ 





cos8(p) sin6(p) 
— sm9(p) cos8(p) 



iprip) 

Tplip) 



being 



tan 26 (p) = 



Uoip) 



(13) 



(14) 



(15) 



Thus the Hamiltonian of the gapped phase, in diagonal form, reads (after dropping an unimportant constant that 
amounts to the ground state energy): 



H LE = &»(P) [= tiip)Mp) ~ 1>t(p)1>v(p) :] 



(16) 



where w(p) = \J ujoip) 2 + m 2 . We associate ipl(p) (iptip)) with the creation operator for particles in the valence 
(conduction) band. 

Before considering quantum quenches, let us briefly discuss some of the the properties of the ground states of the 
Hamiltonians Hq and -Hle- In what follows, these states will be denoted as |$o) an d |3>), respectively. As mentioned 
above, the spectrum of Hq is gapless, and the fermion occupancies in the ground state |<E>o) are: 



n r ip) = (* |$(pM-(p)|*o) = 0{-p), 
mip) = ($o|W(p)^(p)I*o> = 0(p). 



(17) 
(18) 



That is, all single-particle levels with negative momentum are filled out. However, i?LE has a gapped spectrum and, 
therefore, when constructing its ground state, |$), only the levels in the valence band (which have negative energy) 
are filled, whereas the levels in the conduction band remain empty: 



».(P) = W1(P^(P)I*>=1 1 
n c (p) = (*|Vj(p)V'c(p)|*) = 0. 



(19) 
(20) 
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B. Quench from the gapped to the gapless phase 



The first situation we shall consider is when g(t) = g 0(—t) in Eq. (8), so that the spectrum of the Hamiltonian 
abruptly changes from gapped to gapless {i.e. quantum critical). In the following we denote Hi = Hx^{t < 0) = £/le 
and Hf = Hj^{t > 0) = H . The time evolution of ifj r j for t > is thus 

f,i(p,t)=e ±iWo(p) Vr,i(p). (21) 

We first consider the zero temperature quench in this section, and postpone to Section V the discussion of the more 
complicated finite-temperature case. In the zero temperature case the initial state p\_ = |$)($|. Notice that, in this 
state, (3>| cos20(x)|$) = (<&| cos2<£(x)|$) = Rc($\e~ 2lv ^\$) = —(ipl(x)ipi(x)) ^ (the minus sign stems from the 
eigenvalue of the operator T = T) r r}{), whereas in the ground state of Hq, |3?o)j the expectation value of the same 
operator vanishes. Therefore, it behaves like an order parameter in equilibrium, and we can expect that it exhibits 
interesting dynamics out of equilibrium. Indeed, 

C(t) = (e- 2 ^*'*)) = -lj2(ti(pMl(P,t)) = -L^V^to* sm29(p), (22) 

p p 

where, to evaluate the expectation value we have set T — and therefore (ipl(p)ipi(p)) = — |sin2#(p), as it follows 
from Eqs. (14,19,20). The above expression can be readily evaluated by recalling that sin26>(p) = m/u>(p), which 
yields, in the L — > oo limit, 

„/ \ f + °° dp cos 2ujn(p)t / tyi \ , . 1 frn 9mt . „. 

C(t) = m -f- -^-=(- — K (2mt) ~ —\\—e \ 23 

w J 2tt io(p) \2nvJ v ; Av\ id v ' 

where Kq is the modified Bcsscl function. Thus we see that the 'order parameter' (cos2</?(0, t)) decays exponentially 
at long times at T = 0. The decay rate is proportional to the gap between the ground state (the initial state) and the 
first excited state of the initial Hamiltonian Hi = H. The existence of this gap means, in particular, that correlations 
in the initial state between degrees of freedom of the system are exponentially suppressed beyond a distance of the 
order £ c w v/m. Upon quenching, the evolution of the system is dictated by a critical Hamiltonian, Hq, that is, a 
Hamiltonian describing excitations that propagate ballistically along 'light cones' corresponding to the 'trajectories' 
x ± vt. Thus, as discussed by Calabrese and Cardy [9, 15] the correlation length scale characterizing the initial state 
translates into an exponential decay in time of the order parameter at long times. This exponential decay is also 
found (for the same type of quench) in the semiclassical limit of the sine-Gordon model (see Sect. IV below). 
Next we shall consider the (equal-time) two-point correlation function of the same operator, namely: 

g{x,t) = ( e -2^(x,t) e 2^(o,t)^ ( 24 ) 
When using the fermionic representation of e 2 ^^' 1 ^, and upon expanding in Fourier modes we arrive at: 

Q( x ,t) = -^ e^-P^e-* ^ 1 ^ ^ 2 ^ (25) 

Pl,P2,P3,P4 

Applying Wick's theorem, there are three different contractions of the above four fermion expectation value, which 
can be evaluated using Eqs. (14,19,20). This yields the following contractions: 

(4(p)^i(p)) = (Mp)Mp)) =0 (26) 

(^( P )Mp)) = = -isin2%), (27) 

(4(P)MP)) =sin 2 0(p), (28) 

(MP)4(P)) =cos 2 0(p). (29) 
Hence, in the thermodynamic limit (L —¥ oo), we obtain: 

G(x,t) =C(x,t)C(0,t) + \^S(x)+H(x,t)\ 2 (30) 



where 

f + oo 



H(x, t)= JL S mpx^f(e- aoP = —Ki — ^ 

Jo 27T Lo(p) 2ttv \ V J 



(31) 
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in the limit ao — >• 0. Therefore, for \x\ 7^ 0, 

G(x,t) -- 



[[K Q (2mt)} 2 



Ki 



1 2 N 



(32) 



Let us examine the behavior of this correlation function in the asymptotic limit where \x\ £ c = v/m and 2vt ^> £ c . 
Since the Bcsscl functions decay exponentially for large values of their arguments, the leading term in Q(x, t) depends 
on whether t < \x\/2v or t > \x\/2v. Thus, 



G(x,t) 



167TVX 



32nv 2 t 



-2m\x\/v t > \ x \/2v, 
t < \x\/2v. 



-4mt 



(33) 



These results are in also agreement with those obtained using a mapping to a boundary conformal field theory by 
Calabresc and Cardy for general quantum quenches from a non-critical into a critical state. [9, 15] 



C. Quench from the gapless to gapped phase 

We next consider the reversed situation to the one discussed in the previous subsection. In this case, we set 
g(t) = g 9(t) in Eq. (8), i.e. the initial state is critical and corresponds to the ground state of H\ = H , whereas 
the time evolution is performed according to Hi = Hue,. We shall consider the same correlation functions as in the 
previous subsection and therefore it is convenient in this case to obtain the time evolution of the operators ip r (p) and 
ipl(p), whose action on the initial state is known [cf. e.g. Eqs. (17,18)]. Once again, we first restrict ourselves to the 
T = case and defer the discussion of finite temperature effects to Sect. V. We first note that the time-evolved Fermi 
operators can be related to the operators at t = by means of the following (time-dependent) transformation: 

1> r (p, t) = e iH ^ r (P>~ iHtt/H = f(p, t)Mp) + 9*(P, VMP), (34) 
MP,*) = e im Uvy- %mlh - 9*(p,t)A(p) + f*(j>,t)Mp), (35) 

where f(p, t) = cosuj(p)t — i cos 29(p) smuj(p)t and g(p, t) = i sin 28(p) shiuj(p)t. This transformation can be shown to 
respects the anti-commutation relations characteristic of Fermi statistics, and it is therefore a canonical transformation. 
Using Eqs. (34) to (35), we can now compute the decay of the order parameter operator e 2lip ^ x,t \ The calculation 
yields: 

C(x,t) = -| Re LP (P. *)S(P. *)] • (36) 

In deriving the above expression we have used that f(—p,t) = f*(p,t) and g(—p,t) = g(p,t), which fol- 
lows from cos26>(— p) = — cos29(p) because cos29(p) = uj (p)/oj(p). Thus, setting Re[f*(p,t)g(p,t)] = 
— cos 29{p) sin 29(p) sin 2 uj(p)t and taking L — s- +00, we find: 

C(x,t) = 2 / + °° ^^M e -^sin 2 u{p)t = A(ma ) + ^-d(2mt), (37) 
2tt [w(p)J 2ttv 



i) 



where ci is the cosine integral function. The first term is a non-universal constant that depends on the short-distance 
cut-off ao introduced above (cf. Eq. 9). For long times this expression can be approximated by 

C(x,t) ~ A(ma ) + —— sin2mt + 0{t~ 2 ). (38) 

Hence we conclude that, when quenched from the critical (gapless) phase into the gapped phase, the order parameter 
exhibits an oscillatory decay towards a (non- universal) constant value, A(mao). 

Using similar methods the (equal-time) two-point correlation function, Q(x, t) = (e^'^'^e 21 ^ ' 4 '), can be also 
evaluated. The resulting expression can be cast in a form identical to Eq.( 30) of subsect. Ill B. In the thermodynamic 
limit, we find that, in the present case, the function H(x,t) takes the form: 

H(x,t) = -^°°^|-l + [l-cos2 W (p)t]^ I |e^ sinpx (39) 
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We have been unable to obtain an closed analytical expression for this function at all times. However, in t — >• +00 
limit, in which case the term in the integrand proportional to cos2w(p)i oscillates very rapidly and therefore yields 
a vanishing contribution, an analytical expression can be obtained. Upon performing the momentum integral, we 
obtain the following result for large |x| (after taking the limit ao — > 0): 

w <*'*->~>*wjiF- (40) 

Hence, we obtain the following asymptotic behavior of the two-point correlation for t — > 00: 

4v 4 



l 2 



lim Q[x,t) = [A{ma a )Y + j—^ 4 6 . (41) 



This result is clearly different from the equilibrium behavior of the same correlation function in the gapped phase, 
where it decays exponentially to a constant. [31, 32] Instead, when the system is quenched from the gaplcss into the 
gapped phase, we find that both the order parameter and the two-point correlation function (Eqs. 38 and 41) decay 
algebraically to constant (non- universal) values. 

IV. THE SEMICLASSICAL LIMIT 
A. Introductory remarks 

A good approximation to the sine-Gordon model (cf. Eq. 3) can be obtained in the limit where K -C 1, which 
corresponds to the K <C 1 limit in the original notation of Eq. (1). In this limit, we can expand the cosine term 
in (3) about one of its minima, e.g. tp = 0. Retaining only the leading quadratic term yields the following quadratic 
Hamiltonian for the boson field f{x): 



H sG ~ H sc = J 

Z7T 



(a^(*)) 2 + K(d x 0(x)) 2 : + kV9 ®£ 2 I dx : p\x) : . 

2na 1 



Within this approximation, the problem of studying a quantum quench in the sine-Gordon model becomes akin to 
the general problem studied in Ref. 26. To see this, let us first expand •pix) in Fourier modes: 

where tJo(q) = v\q\; the 6-operators introduced above obey the standard Heisenberg algebra: 

[b{q),b\q')] =<W, (44) 

commuting otherwise. The first two terms in Eq. (43) are the so-called zero-modes, whose dynamics is only important 
at finite L. In what follows we restrict our attention to the thermodynamic limit [L — > 00) and therefore neglect the 
dynamics of those zero modes. Introducing (43) into (42), the Hamiltonian takes the general form: 

H(t) =^H [oj (q) + m(q, t)} tf(q)b(q) + \ £ hg(q, t) [b(q)b(-q) + tf(q)b\-q)] , (45) 
1 <? 

with the following identifications oJo(q) = v\q\ and g(q,t) = m(q,t) = 2vg(t)n 2 /\q\a,Q K As in the study of the 
Luther-Emery limit, we shall assume that g(t) = g 9(—t), which corresponds to a quench from the gapped to a 
gaplcss phase, [50], or g(t) = g 0(t), which corresponds to a quench for the gapless to the gapped phase. Following 
the procedure described in the Appendix of Ref. 26, the quench dynamics of this Hamiltonian can be solved by the 
following canonical transformation (indeed, the bosonic version of Eqs. 34, 35): 

b(q,t) = f(q,t)b(q)+g*(q,t)tf(-q), (46) 

where 

f(q,t) = cosuj(q)t -ism.u(q)t cosh2f3(q), (47) 

g{q,t) =ismu(q)t smh2/3{q). (48) 
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Introducing m 2 = Agv 2 ^ 2 /^ K 72 , which is the the gap in the frequency spectrum of the gapped phase and setting 
m(q) = g(q) = m /2u;o(g), the parameter /3(g) satisfies: 

m 2 12 

tanh2/3(g) = 7 (49) 

uj 2 (q) + m 2 /2 

and the frequency: 



u(q) = y/u (q) 2 +m 2 . (50) 
is the dispersion of the excitations in the gapped phase. 



B. Quench from the gapped to the gapless phase 



Let us begin by discussing the situation where g(t) = g 9(—t). In this case, the initial state is the ground state of 
the following Hamiltonian (we omit the zero-mode part henceforth): 



Hi = H sc = ^huj(q) at(q)a(q), 

9# 



(51) 



where the operators a(q) and al(q) are bosonic operators related to b(q) and b'(q) by means of the following canonical 
transformation: 



a(q) = cosh /3(g) b(q) + sinh/3(g) ¥(—q), 



(52) 



with /3(g) satisfying Eq. (49). At t = the Hamiltonian abruptly changes to Hf = Ho, which is diagonal in the 6(g) 
and b'(q) basis, namely, 



H Q =J2^o(q) b*(q)b(q). 

qjtO 



(53) 



In this case the evolution of the expectation value of the order parameter operator e - 2%< t>( x ) = e - lK v>{x) or its 
correlation functions can be obtained from the knowledge of the two-point (equal time) correlation function out of 
equilibrium for the boson field <p{x), i.e. J-(x,t) = (ip(x,t)ip(0,t)) — (ip 2 (0,t)), where the expectation value is taken 
over the ground state of H sc but the time evolution is dictated by Hq. To compute this object, we first insert into 
the expectation value the Fourier expansion of tp{x), Eq. (43) and use Eq. (46). Thus, we arrive at: 



2ttv 



sinh2£(g) cos (qx - 2u (q)t) - e lqx sinh 2 /3(g) - cosh 2 /3(g) 



(54) 



Using this result, let us consider the behavior of the order parameter following the quench. Taking into account that 



e \ ( ^)) : W e see that (ip 2 (0,t)) must be evaluated in closed form using Eq. (54) 



Before performing any manipulation of this expression, it is convenient to subtract the constant (<p 2 (0,0)), which is 
formally infinite (i.e. it depends on the short distance cut-off, eto). Thus, taking the thermodynamic limit where 
L — > oo, we obtain: 



(^ 2 (0,t))-(^ 2 (0,0)) = i 



d(qy) 

uj(q) 



uo(q) 



sin 2 ojo(q)t. 



Inserting the expressions for to(p) and uio(p) in the above equation, we obtain: 

(^ 2 (0,t))-(^(0,0)) = -/(2mt/n), 

where f(z) is defined as: 



3/2 
1 1/2 



(55) 



(56) 



(57) 
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being the Meijer G function. [39] Using the asymptotic expansion for this function, f(z) sa 1 — ^jp, and hence the 
long-time behavior of the order parameter is: 

/ e -2i^(0,t)\ = i e -2iK<p{Q,t)\ _ / e 2i^(0,0)\ e ^(l-irmt) > (gg\ 

We next examine the behavior of the two-point correlation function of the same (order-parameter) operator. 

Q{x, t) = ( e 2*0M) e -M*«M)) = e 4-H*,t) y (59) 

where we have defined T(x,t) = ((p(x,t)(p(0,t)) — (Lp 2 (Q,t)). At zero temperature, with the help of Eq. (54), we find 
that 



2 poo 



m f 

F(x,t)-T{x,Q) = -— J 



-{\-cosqx){l-co$2u {q)t), (60) 



N(<?)] 2 o;(g) 



where 



H*,V) = -\ [°° ^-(l-co S qx)e-« a °, (61) 
1 Jo ^(Q) 

being ao is the short-distance cut-off. Evaluating the integrals: 

•w , . . ,. . , , f[m(x/v + 2t)} + f[m(x/v - 2t)] 

F{x, t) - F(x, 0) = f(mx/v) + f(2mt) - n 2 J 1 w ^ (62) 

where f(z) has been defined in Eq. (57). Thus, asymptotically (for max{|x|/2?j, <} 3> to" 1 ): 



\ x U 2v ,foi t > \x\/2v 
, for t < \x\/2v, 



G(x, t) = e* G(x, 0) x _ kW ^™ (63) 



where Q(x,0) describes the correlations in the initial (gapped ground) state, and exhibits the following asymptotic 
behavior: 

G(x, 0) ~ B(ao) ( 1 - K 2 —e- m W/ v ) (64) 
V rax / 

where B(ao) is a non-universal prefactor. Thus we see that the asymptotic form of Q{x 1 t) (Eq. 63), as well that of 
the order parameter, Eq. (58), have the same form as the results found limit the Luther-Emery limit, and also agree 
with the results of Calabrese and Cardy based on a mapping to boundary conformal field theory. [9, 15] 



C. Quench from the gapless to the gapped phase 



In this case, the system finds itself initially in the ground state of Hi = Ho , and suddenly (at t = 0) the Hamiltonian is 
changed to H{ = H sc . For this situation convenient to obtain the evolution of the observables from the time-dependent 
canonical transformation of Eq. (46), where (3{q) and w(q) are given by Eq. (49) and Eq. (50), respectively. In this 
case, 



(<p(x,t)tp{0,t)) 



1 



E 

9#0 



2ttv \ 



sinh2/3(g) cos (qx ~ 2u{q)t) + e lqx sinh 2 /3(q) + e~ iqx cosh 2 (3(q) 



(65) 



As in the previous subsection, (e 



-2id>(x) 



^ 2 (0,t))-(^ 2 (0,0)) = i 



e \ (v (°>*)> ) and using Eq. (65), we find that: 

2 



d{vq) 



^o(g) 
w(q) 



- 1 



sin u)(q)t 



(66) 



Note, interestingly, that this result can be obtained from Eq. (55) by exchanging u>o(q) and ui(q). However, when 
evaluating the integral we find that (<^ 2 (0, t)) = +oo, for all t > 0, due to the presence of infrared divergences that are 

not cured by the existence of a gap in the spectrum of Hf = H sc . Thus, we conclude that (e" 2 *^)) = e -2 ^"^ 
vanishes at alH > 0. 
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The above result for the evolution of the order parameter seems to indicate that the system apparently remains 
critical after the quench. This is conclusion is also supported by the behavior of the two-point correlation function of 
the operator e 2i ^ x ^: Let Q{x,t) = ( e *<K*,t) e -*<Ko,t)) = e ^r(xd)^ where jr( x ^) = (<p(x,t)<p{Q,t)) - (ip 2 (0,t)}. Using 
Eq. (46) and Eq. (43), we arrive at the following result (at zero temperature, and for L — » +oo): 

F{x,t) - T{x,Q) = ^ [°° d [ vq ) (1 - cosgx) (1 - cos2a,(g)t) (67) 

where 

F{x,0) = - r ^r(l - cosqxy-^ (68) 
Jo Zuj o{q) 

To illustrate the above point about the apparent "criticality" of the (asymptotic) non-equilibrium state, we can analyze 
the behavior of the two-point correlation function, G(x, t), in two limiting cases, for t = and t — > +oo. At t = 0, the 
correlation function, as obtained from Eq. (68), reads: 

g(x,0)=A(a ) (^) 2K , (69) 

where A(ag) depends on the short-distance cut-off a^. Thus, the correlations are power- law because the initial state 
is critical. In the limit where t — > +00, the part of the integral in Eq. (67) containing the term cos 2uj(q)t oscillates 
very rapidly and upon integration averages to zero. The remaining integral can be done with the help of tables, [39] 
yielding: 



F(x, t -> 00) - F(x, 0) = -^Gf 



1 1 
110 1/2 



-04 ^ Av2 

where G 2 , 2 is a Meijer function Using the asymptotic behavior of the Meijer function, [39] we obtain: 



(70) 



lim g(x,t) =B(a ) (— ) , (71) 
\mx ) 

with -B(ao) being a non-universal constant. Thus, although initially the system is critical and therefore correlations 
at equilibrium decay as a power law with exponent 2 k 2 , when the system is quenched into a gapped phase (where 
equilibrium correlations exhibit an exponential decay characterized by a correlation length £ c w v/m), the correlations 
remain power- law, within the scmiclassical approximation. The exponent turns out to be smaller, equal to k 2 , which 
is half the exponent in the initial (gapless) state. In other words, within this approximation, it seems that the 
system keeps memory of its initial state, and behaves as if it was critical also after the quench. This behavior seems 
somewhat different from the results obtained for the same type of quench in the Luther-Emery limit, where both the 
order parameter and the correlations for t + 00 approach a constant value, A(mao) (unless the non-universal amplitude 
A(mao) = 0, which seems to require some fine-tuning). Whether the differences found here between the Luther-Emery 
and the semi-classical limits are due to a break-down of the quasi-classical approximation, which neglects the existence 
of solitons and anti-solitons in the spectrum of the sine-Gordon model, or to a qualitative change in the dynamics as 
one moves away from the Luther-Emery limit, it is not clear at the moment. To clarify this issue will require further 
investigation with more sophisticated methods. 



V. DYNAMICS IN THE LUTHER-EMERY LIMIT AT FINITE TEMPERATURES 



In this section we shall consider that the initial state of the system corresponds to thermal mixed state, which 
describes a sGM system in contact with an energy reservoir (i.e. the canonical ensemble) at a temperature T = 
The state is thus mathematically described by a density operator p\ = e~ Hi / T /Z\. We shall assume that the coupling 
to the thermal bath is turned off at t = 0, and the system subsequently evolves unitarily in isolation according to Hi. 
We shall consider only the Luther-Emery limit of the sGM, where exact results can be obtained at all temperatures 
within the sGM model. The latter is not true in the semiclassical limit discussed above because this approximation 
only captures the breather part of the spectrum and not the solitonic part. We shall therefore not consider finite 
temperature quenches in this limit here. 
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A. Quench from massive to massless 



Consider first the quench from the gapped to the gapless phase. The initial (gapped) Hamiltonian Hi = H^ F is 
thus diagonal in the valence and conduction fermion basis and therefore immediately follows 



(^(pM>(pO) = &[-«(p)]V> 
(ipi(p)i> c (p')} = Mw(p)}6 pp >, 

(^G#c(p')> = o 



Here f F is the Fermi factor and 



= Tr{ 



,-0H< 



}■ 



(72) 
(73) 
(74) 

(75) 



The hnal Hamiltonian Hf 
to (29): 



Hq; thus, using Eqs. (72)-(74) we can compute the Unite temperature versions of Eqs. (27) 



(ip R (p)?p R (p)) = cos 2 8 p f F [w(p)] + sin 2 6 p f F [-uj(p)] 



sm 



J F [u(p)} + cos 2 9 p f F [-uj(p)} 



(^ r {p)Mp)) = (^l(p)Mp)) = -^sin2^tanh 



/3w(p) 



(76) 
(77) 

(78) 



Note that these averages automatically vanish if the fermion operators are evaluated at different values of p because 
of momentum conservation. The time evolution is again dictated by Hf = Hq as in the zero-temperature case. 
Let us next consider some interesting observables. We begin with the order parameter, which reads: 



dp cos 2cl>q (p)t h/3uj(p) 
tcLnli ~~ , 



2tt 



>(P) 



(79) 



This integral can be transformed into an infinite sum by expanding tanh/?e/2 in powers of e ^ e and integrating term 
by term. We thus obtain the following low temperature expansion: 



(80) 



Since the function Kq decays exponentially for large value of its argument, from this expression we see that when 
the temperature is decreased, less terms are needed to approximate the sum. In particular, at zero temperature 
(/3 — > oo) only the n = term contributes and we recover the zero temperature result of Eq. (23). At finite but low 
temperatures, the asymptotic behavior at long times is an exponential decay where the characteristic time decay is 
fixed by the gap. However, at higher temperatures, more terms contribute to the sum whereas the alternating sign 
leads to some partial cancelations. As a result, we expect a faster decay in time of the order parameter. To further 
analyze the long time behavior in this regime, we shall use an identity of Bcssel functions [see Eq. (Al)] which allows 
us to obtain the following high temperature expansion: 



C(t;0) 



2^^ 



cxp 



7T 2 (2/ + 1) 



(h/3m) 2 + tt 2 (2Z + l) 2 



(81) 



Note that, for this expansion, the higher the temperature the smaller the number terms that needs to be retained 
to accurately approximate the sum. In particular, the infinite temperature limit, f3mh <C 1, only the I = terms 
contributes and thus the decay in time is exponential, 



C(t; P) « — exp [-t/r c ] , h/3m <C 1 

ZTTV 



(82) 



but now, the characteristic decay time r c is now fixed by the inverse temperature: 
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Thus, to summarize, the asymptotic behavior of the order parameter following a quench from the gap to the gapless 
phase is described by an exponential decay both at very low and very high temperatures. At very low temperatures, 
the characteristic decay time is given by the frequency gap, m, but as temperature of the initial state is increased, 
the characteristic time is determined by the temperature. The behavior for intermediate temperatures is a crossover 
between these two exponentially decaying behaviors. Moreover, it is also worth mentioning that the exponential 
decay at large temperatures is characteristic of a critical theory at finite temperatures. Indeed, in the sGM we expect 
that, as the temperature is raised well above the gap energy scale, Hm, the the properties of the system will become 
indistinguishable from those of a critical system. 

As for the the finite-temperature correlation function, it can be again cast in the same form as the zero temperature 
case, Eq. (30), with the function t; (3) having the following t — > +oo limit: 



lim W.(x,t;(3) 

t— >+oo 



dp oj (p) 
2tt lo{p) 



tanh 



smpx 



(84) 



Using the same technique as for the order parameter, we evaluate this function as 



K t m 



(85) 



The reasoning is the same as before: at very low temperatures, only the term with n = contributes, and we recover 
the zero temperature expression. However as the temperature increases more terms become important. We can use 
the dual expression Eq. (A2) which applied to Eq. (85) yields 



lim %(x,t\fS) 

t— J- + 0O 



2^E eXP 

H zez 



vnp 



(86) 



This function decays exponentially for long distances (only the I = term contributes at high temperatures), with a 
characteristic correlation length that is fixed by the temperature of the initial state: 



lim H(x,t;P) 

t—t-\-oo 



1 



,-|*|/€ 



2vhf3 



_ vh/3 

7T 



(87) 



Upon introducing this results into (the finite temperature equivalent of) Eq. (30), we arrive at an expression whose 
asymptotic behavior again depends on whether |a;| > 2vt or \x\ < 2vt (with correlation time/lengths given by m or 
T depending on the temperature range). Thus, the correlations at finite temperature also will exhibit the so-called 
'light-cone' effect. [9, 15] 



B. Quench from the gapless to the gapped phase 

We next consider that the system is quenched from the gapless phase intro the gapped phase. Thus we need to 
assume that the system was initially described by H\ = H^e and in contact with an energy reservoir at temperature 
T = /3 _1 . The following expectation values (understood over the initial thermal ensemble) will be required in the 
calculations to follow: 

(^(P)fe(p))=/FN(P)], (88) 
(4(P)^(p))=yt[-Wo(p)], (89) 
(^ R (P)MP)) = (<il>l(p)i> R {p)) = (90) 



and hence 



(ti( P )Mp)) = \( 1 + cos 20 p tanh ^y^) ( 91 ) 

(ri(p)Mp)) = \ (l- cos 28 p tanh (92) 
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Thus, the time evolution of the order parameter can be obtained and reads: 



C(t;(3) 



dp mu (p) h(3uo(p) - 

. .„ tanh e 

2tt u(p) 2 2 



sin 2 u)(p)t 



(93) 



As t — > +oo this function approaches a non-universal constant that depends on the energy cut-off ag and the temper- 
ature. At high temperatures, 



C(t,/3) = A(ma ,(3) + hf3 



lQirvt 3 /' 



■ sin (2mt + f ) 



(94) 



Thus, after the sudden quench at t = from a high temperature state in the critical regime into the gapped phase, 
the order parameter shows an oscillatory decay towards a constant value. However, the exponent of the decaying law 
is different from the decaying exponent in the case of a quench from a zero (or low) temperature state. The whole 
picture is the following: the order parameter exhibits an oscillatory decaying low as t -3 / 2 for times smaller than a 
time scale fixed by the temperature, where there is a crossover to a t -1 / 2 behavior characteristic of low temperatures. 

Concerning the two-point correlation function, it can be again recast as in Eq. (30) with the function H(x,t;/3) 
being given by: 



f°° dp f TO 2 1 

H(x, t;/3) = —l — smpx < — 1 + [1 — cos 2oj(p)t] > tanh 

Jo 2tt I uj{p) z J 



h/3uj {p) 



-a p 



(95) 



At long times, the cosine within the integral oscillates very rapidly yielding a vanishing contribution. The remain- 
ing integral can be evaluated using the Cauchy theorem resulting in an infinite sum over positive odd Matsubara 
frequencies. The sum can be performed, yielding 



H(x, t — ► oo) 



771 1 



— TTX / Vh 



Z Z7T Z ZTT 



1 



2h/3v 



coscch 



/ ttx \ m 



e mx/v tan ^ W C96 1 ) 



where &(x, y, z) is the Lerch function. [ 
frequency term: 

T-L(x, t —> oo) a 



Likewise, the long distance behavior dominated by the lowest (Matsubara) 

Hf3m 



1 

"vhj3 



(h(3m) 2 - 7T 2 



-TTVx/h0 



777 , 

c~ rax 
4t? 



tan ■ 



(97) 



At long distances, this function decays exponentially, but two competing length scales appear: v/m and vH/3/tt. The 
largest sets the characteristic length of the decay. 



VI. LONG-TIME DYNAMICS AND THE GENERALIZED GIBBS ENSEMBLE 



It was recently pointed out by Rigol and coworkers [14] that, at least for certain observables like the momentum 
distribution or the density, the asymptotic (long-time) behavior of an integrable system following a quantum quench 
can described by adopting the maximum entropy (also called 'subjective') approach to Statistical Mechanics, pioneered 
by Jaynes. [40, 41] Within this approach, the equilibrium state of a system is described by a density matrix that 
extrcmizes the von-Neumann entropy, S = — Trplnp, subject to all possible constraints provided by the integrals of 
motion of the Hamiltonian of the system. In the case of an integrable system, if {I m } is a set of certain (but not all of 
the possible) independent integrals of motion of the system, this procedure leads to a 'generalized' Gibbs ensemble, 
described by the following density matrix: 

p g G = ^-e-£^S (98) 

Z gG 

where Z g Q = Tr e _ ^ m Am/m . The values of the Lagrange multipliers, A rn , must be determined from the condition 
that 



(I m ) s G=Tr[p I m ] = {I m ). (99) 



14 



where po describes the initial state of the system, and (• • ■ ) g c stands for the average taken over the generalized Gibbs 
ensemble, Eq. (98). Although pi = |$(t = 0))($(t = 0)| in the case of a pure state, as was first used in Ref. 14, 
nothing prevent us from taking po to be an arbitrary mixed state and in particular a thermal state characterized by 
an absolute temperature, T. In this case, the Lagrange multipliers will depend on T or any other parameter that 
defines the initial state. 

Rigol and coworkers numerically tested the above conjecture by studying the quench dynamics of a ID lattice gas 
of hard-core bosons (see Ref. 14 for more details). One of us showed analytically [11] that correlations of Luttinger 
model also relax to averages taken over this ensemble. This result for the Luttinger model was extended to include 
finite temperature fluctuations in the initial state in Ref. 26. The question that naturally arises then is whether 
the family of integrable models studied in this work (see Eqs. (45), and their fermionic equivalences of Eq. (12) and 
(13)) relax in agreement with the mentioned conjecture. In other words, does the average (0)(t) at long times relax 
to the value (O)gG = Tr p g Q O, for any of the correlation functions considered previously? In what follows we shall 
address this question by analyzing quantum quenches in the sGM at zero temperature. The generalization at finite 
temperatures should be straightforward, as discussed in Ref. 26 



A. Quench from the gapped to the gapless phase in the Luther-Emery limit 



In this case, the type evolution of the system is performed by Hq (cf. Eq. 11), which is diagonal in the operators 
n a (p) = : iptiP^aip) '■ (a = r, I). Thus, the generalized Gibbs ensemble is defined by the following set of integrals of 
motion I m — > I a (p) = n a (p). We see immediately that the fact that this ensemble is diagonal in n„(p) means that 
the order parameter {e~ 2%v ^) s q = (ip r (x)ipi( x ))gG = 0, which agrees with the t —> +oo limit of the order parameter, 
was shown in Sect III to exhibit an exponential decay to zero. However, the two-point correlator of e 2lip ^ has a 
non-vanishing limit for t — > +oo. Thus, our main concern here will be the calculation of the correlation function: 



(e -2i<K*) e 2Mo ))gG = ( V ,t (:r)lM ^t (0)Vv(0)) 



gG 



e i(pi-P2)x 

2^ J2 (^l(Pl)4'l{P2W(p3)i>r(p4))gG 



P1.P2-P3-P1 



(100) 
(101) 



Since the ensemble is diagonal in the chirality index, a, as well as momentum, p, we evaluation of the above expression 
can be carried out by noting that: 



gG 



Tr e - £ pV K> (P)r a > (Ptyt (p)^ Q (p) 



ViP 1 



Tr e 



■5w A a ,(p)/ Q ,(p) 



e A«(p) + I 



(102) 



where the Lagrange multipliers X(q) can be related to the values of the same expectation values in the initial states 
by imposing their conservation, that is, 



(V4(p)Vv(p))gG 



1 



e Mp) + i 
1 



»A r (p) 



= (^( P )Mp)) = cos 2 6(p) 
= <Vj(f#r(p)) =sin 2 0(p). 



Hence, 



(V'J(Pl)V'i(P2)V'; t (P3)V'r(P4)>gG = (V>J(Pl)Vv(P4))gG X (<lpi (p 2 )lpj (P3))gG 

= fipi,pjp2,p 3 sin2 0(Pi) (l ~ cos2 0(P2)) 



U Pl:Pi U P2:P3 



sm 2 6( Pl )sm 2 6(p 2 ). 



Introducing the last expression into Eq. (101) yields: 



(e- 2 ^^e 2 ^ >) gG = 



,tpx 



sin 2 0{p) 



(103) 
(104) 

(105) 
(106) 
(107) 

(108) 



and using that sin 2 6{p) = (1 — cos 2 29(p))/2 and cos2^(p) = u)o(p)/ V^oiP) + m2 ' we nn( ^ (^ or x ^ ®)-> 

■ 



/ m y 
\2ttv) 



(109) 



which is equal to the t —> +00 limit of Eq. (32). 
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B. Quench from the gap less to the gapped phase in the Luther-Emery limit 

In this case the initial state is the gapless ground state of Hq, Eq. (11), whereas the Hamiltonian that performs the 
time evolution has a gap in the spectrum and it is diagonal in the basis of the i/) v (p) and ipc(p) Fermi operators (cf. 
Eq. 16). Therefore, the conserved quantities are 

I v (p) = n v (p)=i)t(j>)i> v (p), (110) 

/c(p) = n c (p)=^(p)Vc(p) (HI) 

The associated Lagrange multipliers (at zero temperature), \ v (p) and X c (p) can be obtained upon equating 
(Iv,c(p))sG = (tf(0)|J„, e (p)|tt(0)>. This yields: 

(I„(p))gG = eX J ) + 1 = *("P) sm 2 B{p) + m cos 2 9{p) (112) 

(Ic(p))gG = eK( l) + i = cos2 W + tf (P) sin2 W > ( 113 ) 

where denotes the step function. Using these expressions we next proceed to compute the expectation values of 
the following observables: 

1. Order parameter 

We start by computing the order parameter, 

(e -*V(x) >gG = = ^]Tsin20( P ) [(J c (p)) gG - (4(p)) gG ] , (114) 

v 

and upon using Eqs. (112) and (113), 

(e - 2M *) )gG = * £8*2*00 0082*00 = - r ^ =Me-^ = A(ma ) 
L ^ Jo 2tt W(p) 

where we have used that cos2#_ p = — cos2# p ; A(mao) is the non-universal constant introduced in Sect. IV C. This 
result agrees with the one obtained in Sect. IV C for the order parameter in the limit t — > +oo. 

2. Two-point correlation function 
We next consider the two-point correlator of the order parameter, namely 

(e 2M*) e 2M0) )gG = 2_ J2 e'^-^(4( Pl )MP2)^(P3)A(P4)) gG . (H6) 

Pl,P2,P3,P4 

The calculation of the average in this case is a bit more involved, but it can be performed by resorting to a factorization 
akin to Wick's theorem. This is applicable only in the thermodynamic limit, as it neglects terms in which the four 
momenta of the above expectation value are equal. These terms yield contributions of 0(1/ L) compared the others. 
When factorizing as dictated by Wick's theorem, the only non- vanishing terms are: 

(l/jj.(pi)lpl(p2)lpi(ps)A(Pi))sG = SpipJp 2P3 (' l Pi{Pl)A(P4,))gG(^HP3) ,l Pl(P2)) g G 

<^(Pl)^(P2)>gG<^ t (P3)V'r(P4)>gG- (117) 

Upon using 

(4( P )Mp)) s g = l#(p) sin 2 2000 + 0("P) f 1 ~ \ ^ 2 %)) - (H8) 
(^( P )Mp))zG = U(~p) sin 2 20(p) + tf(p) (1-1 sin 2 20(p)) , (119) 



(115) 



(V>/0>)Vv(p))gG = ($(pM(p)>gG = -\ sm2d{p) cos2%) sign(p), (120) 
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the average over the generalized Gibbs ensemble of the four Fermi fields on the right hand-side of Eq. (116) can be 
computed and yields the following expression for the two-point correlation function (up to terms of 0(1/L 2 )): 



(e 2M.) e -2M0) )gG 



(121) 



■j ^sin26>(p)cos2%) + \ V e lpx \${p) sin 2 26{p) + ${-p) ( 1 - \ sin 2 26> p ) 
p>0 p L v 7 - 

The first term in r.h.s. of the above expression is just {e 2lip ^) g Q (e~~ 2lip ^)gG (cf. Eq. 115), whereas the first term in 
the right hand-side can be written as 



8(-p) + 2 si e n b) sin2 26 p 



Z2V 

p>0 



-ipx 



-ipx 



sin pa; 



p>0 



w 2 (p) 

2 



lim rl(x,t) 



t— >+oo 



(122) 



(123) 



that is, it coincides with the £ — > +oo limit of the second term in the right hand-side of Eq. 30 in Sect. IIIC (the 
function H(x,t) is defined in Eq. 39). 



C. Quench from the gapless to the gapped phase in the semi-classical limit 

In this case Hamiltonian performing the time evolution is gapless (H ) and thus diagonal in the 6-operators. Hence, 
the conserved quantities are 

7(g) = b\q)b(q) (124) 
The Lagrange multipliers of the corresponding generalized Gibbs density matrix are fixed from the condition: 

< J fo)>«G = p^y— T = (^(0)1^(3)6(5)1^(0)) (125) 
= sinh 2 /3(g), (126) 

where /3(g) is defined by Eq. (49) Hence, using this result we next proceed to compute the order parameter and the 
two-point correlation function. We first note that the order parameter vanishes in the generalized Gibbs ensemble 
since (e- 2 ^) sG = e - 2i ^ 2 (°)>sG an d ((f> 2 (0)) gG = i^(^ 2 (0)) is divergent in the L -t +oo limit (see below). This 
agrees with the result found in Sect. IV B, where it was found that the order parameter decays exponentially in time. 
Thus, in what follows we shall be concerned with the the two-point correlation function. 



1. Two-point correlation function 

Since (e' 2 ^ e 2 ^) gG = e -^r cgG ( x \ where C gG {x) = (^(x)^(0)) gG -(<^ 2 (0)) g G- In order to obtain this correlator, 
we introduce the Fourier expansion of ip(x) (ignoring the zero- mode part), 



2 ^ \u,(q)L 

into the expectation value, and using (126) to evaluate the averages in the generalized Gibbs ensemble, we find that, 
in the thermodynamic limit, 

(v>(sM0)>gG = / cosgz cosh2/3(g), (128) 

Jo 4 woigJ 



and therefore, 



C° G {x) = ( V (x)<p(0)) sG - (^ 2 (0)) gG = - / ^ cosh 2/?(g) (1 - cosgx) (129) 



C(x,0)-— I , , , 2 (l-cosgx) (130) 



m 2 r +oc d(vq) 
TJo u(q) [u a (q)Y 

where C(x,0) = C(x,t = 0) is defined in Eq. (61). Upon comparing the last result with Eq. (60) in the limit where 
t — > +oo, wc sec they arc identical. 
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D. Quench from the gapless to a gapped phase 

In this case the Hamiltonian that performs the time evolution is gapped, whereas the initial state is gapless. Thus, 
differently from the previous case, the Hamiltonian that performs the evolution is diagonal in the a-operators, and 
therefore, the conserved quantities are I(q) = a'(q)a(q). The corresponding Lagrange (at zero temperature) are fixed 
from the condition: 

= P^y-T = (m\aHq)a(q)\m) (131) 
= sinh 2 ^(g), (132) 

where /3(g) is given by Eq. (49). 

In order to obtain the one and two-point correlation functions of e 2l ^ x ^ = e 2%K<p ( x > we first need to write the <p(x) 
field in terms of the a-operators. Upon using the canonical transformation Eq. (52): 




(133) 



Hence, since (e _2 *^ x )) g G = (e~ tKV ^)gG = e~ ^ ^ ^ sG , and (ip 2 (0)) g Q is logarithmically divergent in the thermo- 
dynamic limit (see expressions below), the find that (e~ lK(p ^) g Q = 0. This result is in agreement with the one found 
in Sect. IV C for the order parameter. 



1. Two-point correlation function 

re 2 gG 

Next we consider the two-point correlation function of the same operator, namely [e~ 2l ^ x ^ e 21 ^^ 1 ) g o = e~^ c& ( x ), 
where C sG {x) = (<p(x)<p(Q)) EG - (<^ 2 (0)) gG . Wc first obtain: 

(<p(xM0)) gG = f cosqx cosh2/3(g). (134) 

Jo Mg) 



(135) 



Hence, 

C^ G (x) = - f ^4 cosh /3(g) (1 - cosgz) 
Jo M<?) 

„/ x ™ 2 f + °° diva) , , , 

= C(x,0) + — y 1- cosga; 136 

4 Jo w (g)[o;(g)] 2 

where C(x, 0) is defined in Eq. (68). The latter result agrees with Eq. (67) in the t — > +oo limit. 



VII. RELEVANCE TO EXPERIMENTS 



As we mentioned in the introduction, ultracold atomic systems are the ideal arena to study the quench dynamics 
of isolated quantum many-body systems. This is because they can be treated, to a large extent, as entirely isolated 
systems. Furthermore, since this work is concerned with the quench dynamics of a specific one dimensional model, 
the quantum sine-Gordon model, it is also worth emphasizing that the properties of these systems are highly tunable 
and, in particular, so is their effective dimensionality. Thus, there are already a number of experimental realizations 
of one-dimensional (ID) systems (see e.g. Rcfs. 42-44 and references therein), and in particular, there are also 
experiments where non-equilibrium dynamics has been probed in onc-dimension, e.g. Refs. 6, 45 Thus, there may 
be a good chance that some of the results obtained above may be relevant to current or future quench experiments 
with ultracold atoms. However, since the sine-Gordon model considered in previous sections is nothing but an 
effective (low-energy) description of certain ID systems, any comparison must be done with great care, as there is no 
fundamental reason why the low-energy effective theory should capture the essentials of the (highly non-equilibrium) 
quench dynamics. This is to be contrasted with the equilibrium dynamics, where renormalization group arguments 
show that the sine-Gordon model is indeed sufficient to describe the (universal) properties of certain ID physical 
systems. There is in fact much evidence, both analytical and numerical, accumulated over the years, of the latter 
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fact. However, we are not in a comparable situation in the case of non-equilibrium dynamics, and thus future studies 
should try to address this question more carefully. 

With the above caveat, let us proceed to mention a few situations where the sGM is applicable, at least as a good 
description of the equilibrium state of a system that can be realized with ultracold atomic systems. There are basically 
two kinds of systems, depending on the interpretation of the order parameter operator, e - 2l <t>( x ) . The fist instance is a 
ID Bose gas moving in a periodic potential, where the the sine-Gordon model is the effective field-theory description 
the Mott insulator to superfluid transition (MI to SF) in ID [31, 46]. In this case, the order parameter field is the 
periodic component of the boson density. A quantum quench from the gapped (gapless) into the gaplcss (gapped) in 
this system can be realized by suddenly turning on (off) the periodic potential applied to the ID gas. The evolution of 
the ID density could be monitored by performing in-situ measurements, and the two-point correlations by measuring, 
at different times, the (instantaneous) structure factor using Bragg spectroscopy. 

In the second instance, the order parameter field is interpreted as the (relative) phase of two [47] (or more[48, 49]) 
ID Bose gases coupled by Josephson coupling of two ID Bose gases. Thus, in this setup a quench experiment [16] from 
the gapped (gapless) into the gapless (gapped) phase would correspond to suddenly switching on (off) the (Josephson) 
tunneling, which can be achieved by controlling the (optical or magnetic trapping) potentials that confine the atoms 
to ID. This should be done with care, ensuring the atoms remain in the ID regime both in the initial and final states, 
that is, that e.g. the potential trapping the atoms transversally is always sufficiently tight. The evolution of the 
relative phase can be monitored by analyzing the interference fringes at different times. 

VIII. CONCLUSIONS 

To sum up, we have investigated the time evolution of two-point correlation functions and the order parameter 
after a quantum quench of the relevant operator term in the sine-Gordon model. We considered two different kinds 
of quenches: a quench from a gapped phase to a gapless phase and viceversa. In addition to the initial pure state, we 
considered an initial mixed state coupled to an energy reservoir at finite temperature, that is suddenly disconnected 
at the same time that the quench is performed. In order to compute correlation functions, we studied two limits in 
which the Hamiltonian renders quadratic in terms of either fermion or boson operators, and the dynamics can be 
solved exactly: the Luther-Emery and the semiclassical limits. In the quench from the gaped to the gapless phase, 
the order parameter decays exponentially for long times to the ground state value, with a time scale fixed by the gap. 
In turns, the correlation function exhibit a light-cone effect: it decays exponentially with time with a time-scale fixed 
by the gap until it relaxes to a value that decreases exponentially with distance with correlation length fixed also by 
the mass. These results are valid for both, the Luther-Emery and the semiclassical limits, and agree with the results 
obtained in Ref. 9, 15 . For an initial state at high temperature, these decays are also exponential (and the light-cone 
effect is preserved), but the characteristic time and length is set by the temperature. In between, a crossover connects 
these two limiting cases. 

A general statement for the long time dynamics when the term that opens the gap is suddenly turned on is elusive, 
since in the semiclassical approximation the long times limit of the correlation function exhibit power-law decay with 
distance, and the order parameter vanishes, being these features characteristic of a critical state. On the other hand, 
the results in the Luther-Emery limit indicate a relaxation to a non-universal constant value for long times, signaling 
an ordered state. The latter behavior at the Luther-Emery limit is also present at finite temperature. Whether these 
differences are due to an artifact introduced by the semiclassical approximation or to a very special behavior that 
occurs at the solvable point needs further clarification. 

We have shown that the long time behavior of correlation functions and the order parameter in the different types 
of quenches can be obtained from the generalized- Gibbs ensemble [14] in which the conservation of a certain set 
of independent integrals of motion is fixed as a constrain for the maximization of the statistical entropy. Finally, 
the relevance of the quantum quench dynamics in the sine-Gordon model for cold atomic gases is discussed. The 
supcrfluid-Mott insulator transition appears as the most appropriate scenario to observe the described effects. 
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Appendix A: Some identities involving the Bessel functions 

In this appendix we prove the identities 

exp 



J2 e^Ko (W^ 2 + n 2 ) = tt ^ 



-\i ^a 2 + (2irl + , 



y/a 2 + (2ttZ + , 



E 1 



,m0 



ATi (aV/x 2 + n 2 j = -^exp -[ly/d 2 + (27rZ + 0) 2 



/' 



where </> £ [0,27r). Using the standard integral representation 

1 



K (az) 



dx- 



wc have 
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v k x +k 2 + n 2 + /i 2 



2tt 



/ (ifcj; / lifcj, / dk z S(k z — n) 



a ik x OL+ik z <f> 



k 2 x + k 2 + k 2 z + fi 2 



Next we employ the Poisson summation technique in the form X^nez ^(^z — n) = X)z 



27Ttlk z 



and thus 



riGZ 



]T e*"*# (av^T^) =±-J2 [ dk * [ dk v [ dk *irv 

" iez J J J x v 

= 27 Ey £k k^2 



ik x a-\-ik z <f>-\-2irilk z 



+ k 2 + H 2 



where we introduced the notation 

r x = a, 

Ty = 0, 

r z = <j) + 2irl. 

Using a standard expression for the three-dimensional bosonic propagator, we finally obtain 



e lm *K (oVm 2 + ™ 2 ) = ttE 



where = |r;|. From this equation the first identity Eq. (Al) immediately follows. 

To prove the second identity Eq. (A2) we proceed by integration on both sides of Eq. (Al): 



f daa^e^Ko (ay/ 



fi 2 +n 2 ) = 



■E- 



O 



=K X (ay//x 2 + n 2 ) 
i 2 \ / 



/ daair 



exp 



-M\/a 2 + ( 2nl + 



y/a 2 + (2nl + 0) 2 
-^^cxp [-^v 7 " 2 + (2ttZ + 0) 2 ] 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 
(A6) 



(A7) 
(A8) 



(A9) 
(A10) 
(All) 



(A12) 



(A13) 

(A14) 
(A15) 



Thus, from Eqs. (A13) and (A15), 



- cxp \-^a 2 + (2ttZ + M + C{4>) 



C{4>) is determined from the behavior of Ki(z) for z ^ 0: 



K^z)-- 
z 



and thus, in the limit a — > 0, on one hand, 



lim > i 
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Vm 2 + n 2 



K 1 (ay/fi 2 + n 2 ) = ^ 
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on the other, 



lim — exp —[i\/ct 2 + (2nl + <fi) 2 = — cxp [— [i \2nl 



7T COSh/i(7T — (f>) 

sinh [i-n 



Therefore, both limits coincide and C((f>) = 0, and Eq. (A16) reduces to the second identity. 
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